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Abstract 

A consistent evolution of the strong coupling constant a s from Mz to the GUT 
scale is presented, involving three-loop running and two-loop decoupling. The two- 
loop transition from the MS- to the DR-scheme is properly taken into account. In 
the second part of the paper, the bottom quark mass in the DR- scheme at the 
electroweak/SUSY scale is evaluated with four- loop accuracy. We find that the 
three-loop effects are comparable to the experimental uncertainty both for a s and 
m b . 

PACS numbers: 11.30.Pb, 12.38.-t, 12.38. Bx, 12.10.Kt 

1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is one of the most studied exten- 
sions of the Standard Model (SM). The observation that the gauge couplings of the strong, 
electromagnetic and weak interaction tend to unify in the MSSM at a high energy scale 
f^GVT — 10 16 GeV [1-5] consistent with Grand Unification Theories (GUT) is probably the 
most compelling hint in favour of supersymmetry (SUSY). Constraints from the Yukawa 
sectors of SUSY-GUT models yield appealing predictions for SM parameters, including 
the top quark mass and the ratio of the bottom quark to the tau lepton masses [6-9]. 

Nevertheless, SUSY can only be an approximate symmetry in nature. Several scenarios 
for the mechanism of SUSY breaking, based on very different physical models, have been 
proposed. A possibility to constrain the type and scale of SUSY breaking is to study, 
with very high precision, the relations between the MSSM parameters evaluated at the 
electroweak and the GUT scales. The extrapolation over many orders of magnitude 
requires high-precision experimental data at the low energy scale. A first set of precision 
measurements is expected from the CERN Large Hadron Collider (LHC) with an accuracy 
at the percent level. A comprehensive high-precision analysis can be performed at the 
International Linear Collider (ILC), for which the estimated experimental accuracy is at 



the per mill level. It is obvious that the same precision must be reached also on the 
theory side in order to match with the data [10]. Running analyses based on full two- 
loop renormalization group equations (RGEs) [11-14] for all parameters and one-loop 
threshold corrections [15] are currently implemented in the public programs ISAJET [16], 
SOFTSUSY [17], SPHENO [18], SuSpect [19]. The agreement between the different codes 
is in general within one percent [20] . A first three-loop running analysis, based, however, 
only on one-loop threshold effects, was carried out in Ref. [21,22]. While the three-loop 
running has little effect on the weakly interacting particle sector of the MSSM, this effect 
can be comparable to that of two loop running for strongly interacting particles, as for 
example the squark mass spectrum [22]. 

In this paper, we elaborate an earlier investigation [23] of the evolution of the strong 
coupling a s in MSSM, based on three-loop RGEs [24] and two-loop threshold corrections. 
On the one hand, the three-loop corrections reduce significantly the dependence on the 
scale at which heavy particles are integrated out and on the other hand, they are essential 
for phenomenological studies, because they are as large as, or greater than, the effects 
induced by the current experimental accuracy of a s (M z ). Our aim is to compute a s at a 
high-energy scale \i ~ 0(/xgut) with three-loop accuracy, starting from the SM parameter 
a s (Mz) as input. Additionally, we compare the predictions obtained within the above 
mentioned approach with those based on the leading-logarithmic (LL) approximation 
suggested in Ref. [10]. 

In the context of 50(10) GUT [9,25] the third-family Yukawa couplings may unify 
for a large ratio of the Higgs vacuum expectation values tan/3 ~ 50. Moreover, in all 
SUSY models with large tan /3, the supersymmetric particle spectrum and the Higgs bo- 
son masses are sensitive to the bottom Yukawa coupling [26,27]. In turn, the relation 
between the running bottom mass and the bottom Yukawa coupling is affected in such 
theories by large supersymmetric radiative corrections [28] . For the calculation of higher 
order corrections in the framework of SUSY theories, the DR scheme, based on regu- 
larization by dimensional reduction (DRED) and modified minimal subtraction, is most 
convenient. On the other hand, the values of the SM parameters are presently extracted 
from the experimental data for the MS renormalization scheme. It is obvious that trans- 
lation formulae from one scheme to the other to the appropriate order in the perturbative 
expansion have to be employed for phenomenological studies. 

Owing to these considerations, a precise determination of the running bottom mass in 
the DR scheme is of great importance for SUSY theories sensitive to the bottom quark 
Yukawa coupling. Ref. [29] proposed an approach to the determination of mf R (M z ), 
which avoids the infrared sensitivity problem related with the use of the bottom-quark 
pole mass. The procedure consists of employing the RGEs of QCD and the relation 
between MS and DR masses in order to calculate m^ R , using as input parameter the 
accurately known //& = mjf s (/i{,). A second purpose of this paper is to extend the study 
of Ref. [29] to three- and even to four-loop order [30,31], in order to take full advantage of 
the recent determination of mj^ s (/i{,) with four-loop accuracy [32]. Furthermore, we study 
the phenomenological significance of evanescent couplings related with the application of 
the DR scheme to QCD. 
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The remainder of the paper is organized as follows: in Section 2 we discuss the evalu- 
ation of the strong coupling at the GUT scale from the knowledge of a s (Mz) and propose 
a consistent prescription based on three-loop running and two-loop decoupling. Section 3 
deals with the running of the bottom quark mass from /i — /if, to a decoupling scale of 
super symmetric particles. 



2 Running and decoupling of ct s 

The value of the strong coupling constant, measured at the mass of the Z boson M z , is a 
central quantity in high energy physics. Usually, what is being quoted in the literature as 
a s (Mz) is its value in a theory with five quark flavours, renormalized in the MS scheme. 

In order to be precise, we will refer to this quantity as af S '^ 5 \M z ). 

On the contrary, what is meant by the value of q; s (/xgut) is the strong coupling constant 
in a supersymmetric theory, renormalized in the DR scheme. We will denote this quantity 
by q;? R '^ u11 ^(//gut) in what follows. It is the purpose of this section to study the relation 

between «f S, ' 5 '(M^) and a? R ^ full ^(//GUT) a t three-loop level. This requires the consistent 
combination of (a) the renormalization group evolution of a s ; (b) the transition from the 
MS to the DR scheme; (c) the transition from five-flavour QCD to the full SUSY theory. 
All these issues will be discussed in detail in what follows. 



2.1 Renormalization group evolution 

The energy dependence of the strong coupling constant is governed by the RGE 

" n>0 

The f3 function depends on the underlying theory and on the renormalization scheme. In 
QCD with rif quark flavours, it is known in the MS scheme through four loops [33,34]; 
the first three coefficients read 

* = i (" - h) ■ 3 ' = s( 102 -f"', 

* = 64 ^ " if n ' + ^ n 'l] 

The QCD (3 function in the DR scheme involves evanescent couplings (see Section 2.3 
below) at three loops and higher. Explicit results through four loops are given in Refs. 
[30,31]. 

In SUSY-QCD, the (3 function has been evaluated in the DR scheme through three 
loops from the NSVZ formula [35,36] by a proper shift in the strong coupling constant [24]. 
For six (s) quark flavours one finds 

3 7 S47 

A = j, A = - g , A = - M . (3) 
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2.2 Decoupling of the SUSY particles 

For mass independent renormalization schemes like MS or DR, the decoupling of heavy 
particles has to be performed explicitely. In practice, this means that intermediate ef- 
fective theories are introduced by integrating out the heavy degrees of freedom. The 
parameters of the various effective theories are related by so-called decoupling (or match- 
ing) constants. 

When going from a high to a low energy scale, one may separately integrate out ev- 
ery particle at its individual threshold. This multi-scale approach (sometimes denoted as 
"step approximation") is suited for SUSY models with a severely split mass spectrum [37]. 
As pointed out long time ago, the intermediate effective theories with "smaller" symme- 
try raise the problem of introducing new couplings, each governed by its own RGE. 
This is the case in the Yukawa sector if the SU(2) symmetry is broken, for example. It 
also happens for the "evanescent" couplings, related with the application of DRED to 
non-supersymmetric theories. A comprehensive multi-scale decoupling procedure in the 
context of the MSSM is currently not available, although it could be of phenomenological 
relevance [38]. 

On the other hand, in SUSY models with roughly degenerate mass spectrum at a scale 
M, one can consider the MSSM as the full theory that is valid from the GUT scale /xgut 
down to M, which we assume to be around 1 TeV. Integrating out all SUSY particles 
at this common scale, one directly obtains the SM as the effective theory, valid at low 
energies. As already mentioned above, the strong coupling constants in SUSY and in QCD 
are then related by a decoupling constant ( s . The transition between the two theories can 
be done at an arbitrary decoupling scale /x: 1 

a ™.(»/)( |i ) = C (»/) a DE,(Wii)( |X ). (4) 

( s depends logarithmically on the scale //, which is why one generally chooses fi ~ M. In 
Eq. (4), rif — 6 means that only the SUSY particles are integrated out, while for rif — 5 at 
the same time the top quark is integrated out. This procedure, also known as "common 
scale approach" [38], has the advantage that it avoids the occurrence of intermediate 
non-supersymmetric effective theories. Most of the present codes computing the SUSY 
spectrum [17-19] follow this approach by applying the one-loop approximation of Eq. (4) 
and setting rif — 5 and /i = Mz- 

In this paper, we apply the corresponding two-loop approximation, as it is required by 
a consistent three-loop evolution of the coupling constant [39]. Furthermore, in contrast 
to most analyses performed up to now, we allow for a general decoupling scale \i = /idee 
where the heavy particles are decoupled. 

The decoupling coefficients of Eq. (4) have been evaluated through two loops in 
Eqs. (26) and (29) of Ref. [23]: 

Ci n/) = 1 + - C ( i /} + (-) 2 &2 f) + 0(al) , (5) 

1 In principle the decoupling constants should also carry a label MS or DR. However, in this paper 
decoupling is always performed in the DR scheme, and therefore we omit the label. 
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where a s = a? R '^ uII ^(//), and 
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M is the common mass of the SUSY particles, m t is the top quark mass, and 
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Although the formula for in Eq. (6) has been obtained in the limit m t M, it has 
been shown that it holds at the per cent level for all values of m t < M [23]. The exact 
limit m t = M is also given in Ref. [23]. 



2.3 MS— DR conversion and evanescent couplings 

In order to use Eq. (4), one needs to transform the input value for a s from the MS to 
the DR scheme. The corresponding relation has been evaluated through three loops in 
Ref. [31]. Here, we only need the two-loop expression [30]: 



a 



.\is 



a 



DR 




+ 



127T 2 



n f + 



(8) 



where a» R = a£ R ' (n/ V) and af s = af S ' (nf) (fi) . 

a e = a^^/i) is one of the so-called evanescent coupling constants that occur when 
DRED is applied to non-supersymmetric theories (QCD in this case). In particular, it 
describes the coupling of the 2e-dimensional components (so-called e-scalars) of the gluon 
to a quark. Other evanescent couplings that occur in QCD are r\ r (r = 1,2,3), related 
to the vertex of four e-scalars (see, e.g., Ref. [31] for details). In the following, we will 
assume that QCD is obtained by integrating out the heavy degrees of freedom (squarks 
and gluinos) from SUSY-QCD. In this case, the evanescent couplings can be related to 
the gauge coupling a s as follows: 2 



a 



(full) 

e 

(full) 



(/i)=a DR,(fun )(/i)=?? (f u ii) (/i)) 



(9) 



2 Assuming that DRED preserves supersymmetry. 
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The evanescent couplings in n/-flavour QCD, i.e. a e and r]r are then obtained by 
decoupling relations analogous to Eq. (4). 



2.4 Decoupling for oc e 

For the one-loop decoupling relation that connects the values of a e in five-flavour QCD 
and in SUSY-QCD, we find by a calculation similar to Ref. [23]: 



(full) 

DR,(full) 



Ce = 1 + 



7Y 



m~ — 3m? 



2 



+ C A 



(m|(l + L 5 )-m|(l 



4 (to- 
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(10) 



where we have used the short-hand notation 

,2 



In 



/r 



mt 



x e {t,g,q}. 
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Furthermore, the mixing angle between the squark mass eigenstates to^i and m^ 2 has been 
set to zero. Note that Eq. (10) depends on a specific (s)quark flavour q. This means that 
actually Ce and therefore also a e should carry an additional label q: 

^ {5) (^ = C«f u) ^). 



For example, a e is the coupling constant for the uue vertex, where the decoupling 
constant Q depends on m„. Therefore, if all squark masses are different, one needs a 
separate evanescent coupling for each quark flavour. We remark that gauge invariance 
prevents the occurrence of such a flavour dependence for a^ R '^ 5 \ If all squarks have the 
same mass m ? -, then all evanescent couplings are the same, a. 



«,(5) 



a. 



d,(5) 



and one can drop the flavour label. In the limit m~ g = rriq = M, Eq. (10) reads 



Ce = 1 + 



a DR,(full) 



7T 



l -L t 
2 



+ C A ( -L u ) + C 



(12) 



a e occurs typically at higher orders, and in fact, as we will now argue, one can easily 
eliminate it perturbatively from the three-loop running of a s . 
Because of Eq. (9), one can write 



™/)( /i )=^/)a? R ' ( ^ ) ( /U ), 



(13) 



with 



u(n/) 
>e 



DR,(n/) / s 

1+ a. W^ l+0(aS 

71 



(14) 



where the dependence on a s ,(u ' has been eliminated by using Eq. (4). In case all heavy 
particles have a common mass M, the quantity ^[ n ^ only depends on ln(/i 2 /M 2 ) and thus 
for ji = /idee ~ M no numerically enhanced perturbative coefficients appear. This is even 
true for a mass spectrum where all mass ratios are within the range of one or two orders 
of magnitude. This absence of artificially large logarithms allows us to write 

^W) = «?* (B 'W) + 0(«3 . (15) 
In this way, we can replace a e by af R in the two- loop term of Eq. (8). 



2.5 Evaluation of ck s (/xgut) from ct s (Mz) 

There are various ways to go from ctf S ^ ai \M z ) to q;? R '^ uII ^(/Xgut)- They most impor- 
tantly differ in the degree to which evanescent couplings are needed. We propose the 
following method: 

a s (M Z ) -> Oi s (//dec) -»■ Ols (/idee) 

(16) 

( - } ctf^W) n c* s DR «(//gut) ■ 

The individual steps require 
(i) (3(a s ) in QCD through three loops [Eqs. (1), (2)] 

(ii) the MS-DR relation through order a 2 s [Eq. (8)]; due to Eq. (15) one can set a e = a^ 15, 
at this order 

(Hi) decoupling through order a 2 s [Eqs. (4) with rif — 5, Eq. (6)]; no evanescent couplings 
appear at this order (for a "common scale approach"). Note that if the decoupling 
of the supersymmetric particles were performed in several steps, a e would appear 
in the two-loop decoupling coefficients 

(iv) (3(a s ) through three loops in SUSY [Eqs. (1), (3)] 

An advantage of this procedure as compared to a multi-scale approach is that the RGEs 
are only one-dimensional (there is only one coupling constant), and that for a e one can 
apply Eq. (9) and Eq. (15). 

Let us remark that in principle it is possible to decouple the top quark separately: 

a™.< B >(M z ) n aT S > {5) (»t) ^ a™' (8) (*) ^ af~ S ' {6) W) (") 

and continue with steps (ii)-(iv) of Eq. (16), but now using rif — 6 in all the formulas. 
fit is an additional decoupling scale to be chosen of the order of m t . The only new 
ingredient needed is the decoupling constant for going from five to six quark flavours 
in the MS scheme in step (ii'). It has been evaluated through three and four loops in 
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Refs. [40] and [41,42], respectively. In any case, for a mass spectrum as given by the 
benchmark point SPSla' [10], for example, the separate decoupling of the top quark 
implies a numerically small effect. This can also be established by comparing "Scenario 
D" and "Scenario C" in Ref. [23]. 

For a direct application of Eqs. (4), (6), and (8) to steps (ii) and (Hi) of Eq. (16), these 
equations need to be inverted. In fact, it may be convenient for practical phenomenological 
analyses to combine steps (ii) and (Hi) of Eq. (16) into a single formula: 



a 



DR,(full) = a MS '( n /) <; i 



+ 



MS,(n f ) , 

4 ^ 



7T 



MS,(n,)' 



IT 



ii _ V± - i A»t) _i_ o ( A"f)\2 _ An,) 
8 12 2 2 



(18) 



where the scale /id ec has been suppressed in the notation. The coefficients and 

are given in Eq. (6). The numerical deviation of Eq. (18) from the two-step procedure is 

well below 0.2%. 



2.6 Numerical results 

The result for a£ R ' (full V G UT = 10 16 GeV), obtained using M z = 91.1876 GeV and 

m t = 170.9 ±1.9 GeV, af E ' (5) (M z ) = 0.1189 , M = m~ q = m- g = 1000 GeV (19) 

as input parameters is shown if Fig. 1. The dotted, dashed and solid line are based on 
Eq. (16), where n-loop running is combined with (n — l)-loop decoupling, as it is required 
for consistency (n = 1,2,3, respectively). We find a nice convergence when going from 
one to three loops, with a very weakly /id ec -dependent result at three-loop order. 

For comparison, we show the result obtained from the formula given in Eq. (21) of 
Ref [10]. It corresponds to the dash-dotted line in Fig. 1. In this case, the evolution is 
done in two steps: 

af s ^(M z ) a^^Mz) ^ o^^/iout) ■ ( 20 ) 

In (i), af S ' ( ' 5 \Mz) is converted into a^ R '^ uir> (M z ) by means of resummed one-loop contri- 
butions originating from both the change of scheme (from MS to DR) and the decoupling 
of heavy particles. In step (ii), the evolution from M z to /xgut is performed using the 
one-loop MSSM RGE in the DR scheme (f3$ from Eq. (3)). A remarkable feature of this 
combination is that the /!d cc ~dependence drops out explicitely. However, the difference 
between our three-loop result with two-loop decoupling (upper solid line) and the one-loop 
formula given in Ref. [10] exceeds the experimental uncertainty by almost a factor of four 
for sensible values of /idee- This uncertainty is indicated by the hatched band, derived 
from 5a s (M z ) = ±0.001 [43]. The formulae of Ref. [10] should therefore be taken only 
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Figure 1: a s (/xcuT) as a function of /id ec - Dotted, dashed and solid line: prescription 
discussed in Eq. (16) for one, two, and three loops, respectively. The hatched band denotes 
the uncertainty from the input value a s (Mz)- Dash-dotted line: one-loop running and 
one-loop decoupling as described in Ref. [10]. 



as rough estimates; once precision studies are required, one should rely on the consistent 
treatment of running, decoupling, and MS-DR conversion, as it is outlined here. 

In Fig. 2 we show o; s (/zgut) as a function of M where //dec = M has been adopted. 
Dotted, dashed and full curve correspond again to the one-, two- and three-loop analysis 
and the uncertainty form a s (Mz) is indicated by the hatched band. One observes a 
variation of 10% as M is varied between 100 GeV and 10 000 GeV. This shows that 
the actual SUSY scale can significantly influence the unification, respectively, the non- 
unification behaviour of the three couplings at the GUT scale. For definite conclusions 
also the influence of M on the electro-magnetic and the weak coupling has to be carefully 
studied, of course. 
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Figure 2: o; s (/zgut) as a function of M where /idee = M has been chosen. The notation is 
adopted from Fig. 1. 



3 Evaluating rab(/x) in the DR scheme 

Large values of tan (3 could enhance the importance of the bottom Yukawa coupling in 
the MSSM enormously as compared to the SM. If predictions within this model are 
calculated in the DR scheme, then it is certainly necessary to know mf R with the highest 
possible precision. Experimental data concerning the measurement of the bottom mass 
are typically converted into //& defined by the recursive equation 

fi b = mf S ' (5) (fi b ) . (21) 

In the following, we will provide a precise relation between this value and n\ R '^\fis), 
where fi$ may be an energy scale between the electroweak gauge boson masses and a few 
TeV. _ 

The quantity n\ K,< ^\Mz) has also been considered in Ref. [29] at two-loop level. We 
will comment on the differences between this and our result at the end of this section. 
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(5) 

3.1 Evolution and scheme conversion of ra), 

We follow two different methods. The first one is given by the chain 

MS,(5)/ \ (*) MS,(5)/ \ (»*) DR,(5)/ \ 

The RGEs of QCD required for step {%) are known to four-loop accuracy [44,45], and 
three-loop corrections are available [30] for the MS-DR conversion relation in step {%%) . A 
peculiar feature of this mass conversion is that the evanescent coupling a e occurs already 
at one-loop level: 



m DR = m Ms/ 1 _^ £ + 

V 3tv 

where the dots denote higher orders in a e , r\ r and a s . 
The second method consists of the following sequence: 



•) , (23) 



-r S ' (5) M - -C' (5) M ^ "f^W (24) 

Again, three-loop terms are known for the relation used in step {%'). The running in the 
DR scheme in step (W) involves the non-diagonal RGEs of a^ R '^ and the evanescent 
couplings and r]i 5 \ The corresponding (3 functions are known to four-, three-, and 
one-loop order, respectively [30]. 

Because of these multi-dimensional RGEs required for Eq. (24), the sequence of 
Eq. (22) is preferable to the latter in practical calculations. However, it may be inter- 
esting to see how the perturbative series behaves for the two methods, and to which 
extent they are numerically equivalent. In particular, this may provide an estimate of the 
uncertainty in the relation between n\ S '^\jj,i)) and m^ R '^(//s). 

Just as for the strong coupling a s , the consistent evolution of requires that n-loop 
running is combined with (n — l)-loop MS-DR conversion. Assuming that the couplings 

af S '' 5 '(M^), a:? R '^(Mz), af\Mz) and rji 5 \Mz) are known, all operations are within 
five-flavour QCD and do not require any decoupling. This allows it to apply Eq. (22) at 
four- loop order. This is not strictly the case for Eq. (24), because the (3 functions for the 
evanescent couplings are not known to four loops. We ignore this fact when evaluating 
the three- and four-loop evolution along Eq. (24), expecting these effects to be small. 



3.2 Input values for a e and rj r 

In addition to mjf S ' ( ' 5 ' ) (/! ;) ), the only input value for the determination of m° R '' 5 Vs) i s 
cxf S '^\Mz). However, either of the two methods of Eqs. (22) and (24) requires to switch 
to the DR parameters a? R,< - 5 ' ) , ati\ and r]^ at some point. 

In principle, any choice of a£ R ^ 5 \M z ), af\M z ) and t]i 5 \M z ) is allowed which obeys 
Eq. (8) for the given value of «f S, ' 5 '(M z ). Any such choice would simply correspond 
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to a particular renormalization scheme. But let us assume that QCD is the low energy 
effective theory of SUSY-QCD, and thus Eq. (9) holds. Then all DR couplings of five- 
flavour standard QCD are uniquely determined by the corresponding decoupling relations. 
In contrast to Section 2, we now want to go to the four-loop level, and thus we take the 
one-loop decoupling relation for a e into account 3 , see Eq. (10). The r) r , on the other 
hand, occur only at four- loop level and can be decoupled trivially by setting 7]r nU \ndec) = 

??r 5) (/idcc)- 

Since a° R '' M1 ' is not known a priori, one cannot use Eq. (10) directly in order to derive 
af\ Rather, we start with a trial value for a? R ^ full ^(/Xdec) an d obtain the corresponding 
o^e\l^dec) through Eq. (10), as well as a? R '^(/x dec ) through Eq. (6). Then we evaluate 
ctf S ^ 5 \fidec) through Eq. (8), and from that af S ' < ~ 5 \M z ). The trial value for a? R '^ full ^(/Xdec) 
is systematically varied until the resulting af S, ' 5 '(Mx) agrees with the experimental input. 

As indicated above, current knowledge of the renormalization group functions allows 
us to follow the sequences Eq. (22) and Eq. (24) at four-loop level. Strictly speaking, this 

requires also that q^ K '^\m z ), a?'(M z ), and r]i 5 \M z ) are derived from af S '^ 5 \M z ) at 
four-loop level. The decoupling of the SUSY particles mentioned above should thus be 
performed at three-loop level, which is currently not known for any of the three cou- 
plings. However, we may consider q^ K, ^\m z ), af\M z ), and r]f > (M z ) as input, and 

their derivation from af S '^\M z ) only as a guideline, allowing us to neglect higher order 
decoupling effects. We do not expect that a fully consistent evaluation of these parameters 
will change our numerical results for mf R significantly. 

3.3 Numerical results 

Taking Eq. (19) as well as m^ S '^(//j,) = 4.164 GeV [32] as input values and following 
steps (i) and (ii) of Eq. (22), we find (the numbers refer to /idee = Us — M z ): 

(i) The difference between two- and three-loop running on tt\ R ' ( ' 5 \m z ) is 23 MeV, 
while going to four loops further modifies the result by less than 2 MeV. 

(ii) The one-loop term in the MS-DR transition amounts for fx = M z to about 29 MeV, 
while the two-loop term is only about 2 MeV. 

Following steps (i') and (if) of Eq. (24), on the other hand, one finds: 

(i') The one-loop term in the MS-DR transition amounts for fx — fx^ to about 66 MeV, 
while the two-loop term is only about 6 MeV. 

(ii') The difference between two- and three-loop running on n\ R '^ 5 \M z ) is 18 MeV, 
while going to four loops further modifies the result by less than 0.5 MeV. 



3 Since we take the three- loop beta function for a e into account we would need the two- loop corrections 
at this point. However, we assume that the numerical effect is small. 
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Figure 3: mf R (M z ) as a function of /id ec f° r two-, three- and four-loop running (dotted, 
dashed and solid lines). The upper (lower) curves correspond to the prescriptions (i) and 
(ii) ((i') and (ii')). The band describes the uncertainty from a s or mj (see text). 



In Fig. 3, mf R (Mz) is shown as a function of fid ec where the dotted, dashed and solid 
lines correspond to two-, three- and four-loop running. 4 At each loop level, the upper curve 
corresponds to the path given in Eq. (22), the lower one to Eq. (24). The smaller hatched 

band reflects the uncertainty induced by the current experimental accuracy of a^ IS '^ 5 ' ) (Mz), 
equal to 5a s = 0.001 [43], and the larger one corresponds to 5mb = 25 MeV [32]. The total 
uncertainty from both and a s is thus obtained by adding in quadrature the individual 
uncertainties. 

As expected, the difference between paths Eq. (22) and Eq. (24) decreases as higher 
order corrections are included. At four-loop level, the lines are practically indistinguish- 
able. One notices a slightly better convergence of the perturbative terms when following 
the path of Eq. (24). Note also that both two-loop curves lie on the upper edge or even 

4 For the one-loop result the difference between and (i'),(W) is small because no non-trivial 

matching terms are included at this order. The resulting value mf R (Mz) « 3.05 GeV lies well outside 
the plot. 
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A^S 


A'-dec 


M z 




[I), 


1 


2 


3 


4 (running) 




0.1189 


0.0968 


0.0995 


0.0995 


KH'oJ 


0.1189 
3.055 


0.1200 
2.859 


0.1202 
2.838 


0.1202 
2.836 


m b \Ps) 


3.055 


2.830 


2.807 


2.805 


(i% (ii>) 




a { e\lJ>b) 


0.1649 


0.1501 


0.1537 


0.1538 


DR,(5)/ \ 


0.2153 


0.2316 


0.2336 


0.2343 




4.164 


4.098 


4.092 


4.094 




3.055 


2.823 


2.805 


2.805 








A^dec 


2M Z 




(i), (ii) 


1 


2 


3 


4 (running) 




0.1133 


0.1006 


0.1014 


0.1014 


DR,(5)/ \ 


0.1189 


0.1200 


0.1202 


0.1202 


-r i(5) (^) 


3.055 


2.859 


2.838 


2.836 




3.055 


2.829 


2.807 


2.804 


(</), («') 






0.1590 


0.1545 


0.1559 


0.1560 




0.2153 


0.2316 


0.2336 


0.2343 


m 6 DR ' (5, (») 


4.164 


4.096 


4.091 


4.092 


m t DR ' (5 V S ) 


3.055 


2.822 


2.804 


2.804 



Table 1: Numerical results for the bottom quark DR mass for = Mz using /id cc = Mz 
and /idee = 2Mz, respectively. For convenience intermediate results for a s and a e are 
given resulting from the one-, two-, three- and four-loop analysis. 



outside the uncertainty band of a s and m&, thus proving the importance of the three-loop 
terms. Detailed numerical results for the couplings and bottom quark masses are given 
in Tab. 1. We consider the choices /is = /-tdec = M z , as well as fi$ — = M Z - 

It is worth mentioning that the numerical effect of the decoupling of a? R ' <n/ ' ) (cf. 

Eq. (4)) and a e (cf. Eq. (10)) on m° R,{5) (M z ) is not negligible. For (i nf) = 1 and ( e = 1, 
we observe a decrease of more than 6 MeV for both Eqs. (22) and (24), employed at two-, 
three-, and four-loop order. 

In Fig. 4 we show m^ R (/i) varying /i between 50 GeV and 5 TeV where the two-, 
three- and four-loop results corresponding to the steps (i) and (ii) of Eq. (22) are plotted. 
The decoupling scale is set to /idcc — M — 1 TeV. In Tab. 2 we show for some selected 
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jlx (GeV) 

Figure 4: mP R (/i) as a function of /x. The notation is adopted from Fig. 3. 



H (GeV) 




91.1876 


2.804(16)(20) 


350 


2.528(17)(18) 


500 


2.467(17)(18) 


800 


2.394(17)(17) 


1000 


2.361(17)(17) 


2000 


2.268(17)(16) 



Table 2: mf for various values of the renormalization scale. The two numbers given 
in the round brackets correspond to the uncertainties induced by <5a^ s = 0.001 (first 
number) and <5m^ IS (/i b ) = 25 MeV (second number), respectively. 



values of /i the corresponding result for mf together with the uncertainties arising from 
feMS = o.ooi and 6m^(ji b ) = 25 MeV. 5 

5 In the analysis of Ref. [32] Sa^ s — 0.002 has been adopted, leading to the uncertainty of 25 MeV for 
m&. Still, for illustration purpose we choose i5a^ s = 0.001 in order to obtain the numbers in Tab. 2. 
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Let us now comment on the earlier two- loop calculation of mf R by Baer et al. [29]. 
Our results represent an improvement with respect to several issues. First of all, we 
increased the accuracy of the result by including the three- and the four-loop terms in 
the derivation of mf R . In particular the three-loop terms turn out to be numerically very 
important, while the four- loop result indicates a nice stabilization of the perturbative 
expansion. In addition, we have made a significant conceptual generalization, related 
to the evanescent couplings. Ref. [29] sets a s (Mz) = a e (Mz) in Eq. (23). As we have 
argued above, this corresponds to a particular renormalization scheme within DR, and the 
value for mf K is associated with this scheme. In contrast to that, we have systematically 
derived a e by assuming that squarks and gluinos decouple at a scale /idee and have studied 
the dependence of mP R on this scale. Finally, we have given explicit numerical results for 
the method described by Eq. (24). Ref. [29] correctly claims that the difference to the 
method of Eq. (22) is small, but no numbers or any details of the calculation are given. 
In particular, it is unclear which value was used for a e (/if,) in Eq. (23). 

It should be pointed out that the formulae and prescriptions of Ref. [29] have been 
partially taken over in the outline of the SPA project [10]. However, as in the case of the 
strong coupling, various orders of perturbation theory have been combined inconsistently. 
Thus, let us stress again that the formulae of Ref. [10] should not be taken over literally 
in phenomenological analyses if precision is of concern. 

4 Conclusions 

Application of DRED to non-SUSY theories is rather cumbersome compared to DREG 
because of the occurrence of evanescent couplings. However, when low energy precision 
data such as a^ s {Mz) or m^ LS (fj,i > ) are related to their counterparts in a SUSY theory at 
high energies, one may need to switch between the MS and the DR scheme at some energy 
scale. The conversion formulae will involve evanescent couplings, sometimes already at 
one-loop level, as in the case of the quark mass. 

We have used recent three- and four-loop results for the (3 functions, the quark 
anomalous dimension, and the decoupling coefficients in order to derive a^ R (/iGUT) and 

m^ R '^(/i) from af s (Mz) and m^ IS (/if 1 ) at three- and four-loop level, respectively. 

It turns out that the three-loop terms are numerically significant both for a s and for 
nib. The dependence on where the SUSY spectrum is decoupled becomes particularly flat 
in this case. The theoretical uncertainty is expected to be negligible w.r.t. the uncertainty 
induced by the experimental input values. 

Comparing our results and methods to the literature, we find that the issue of evanes- 
cent couplings has either been ignored (by assuming a e = a s ) or circumvented by decou- 
pling the SUSY spectrum at fi^cc = M Z - We find that at one- and two-loop level, this 
choice does not allow for a good approximation of the higher order effects, if one assumes 
the SUSY partner masses to be of the order of 1 TeV. 

To conclude, we strongly suggest that phenomenological studies concerning the impli- 
cations of low energy data on Grand Unification should be done at three-loop level, and 
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that decoupling effects and evanescent couplings are properly taken into account. 
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